Stepwise quantum decay of self-localized solitons 
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The two-phonon decay of self-localized soliton in a one-dimensional monatomic an- 
harmonic lattice caused by cubic anharmonicity is considered. It is shown that the decay 
takes place with emission of phonon bursts. The average rate of emission of phonons 
is of the order of vibrational quantum per vibrational period. Characteristic time of 
the relaxation is determined by the quantum anharmonicity parameter; this time may 
vary from a few (quantum lattices, large anharmonicity) to thousands (ordinary lattices, 
small anharmonicity) of vibrational periods. 

PACS numbers: 63.20Kr, 63.20Pw, 63.20Ry 

The dynamics of strong nonlinear excitations in polymers and quasi-one-dimensional biomolecular 
chains is an active research field. Theoretical studies of anharmonic perfect lattices have shown the 
existence of localized vibrations (self-localized solitons, SLSs) with frequencies above the phonon band 
or in the gap of the phonon spectrum (seetrEj and references therein) . SLSs are soliton-like excitations 
in discrete lattices, and are thus closely related to ordinary solitons. The existence of stationary and 
moving SLSs was derived in the frame of classical mechanics. Until now there was a little discussion 
about the influence o£ quantum (and thermal) fluctuations on stability of the SLS. An exception is the 
paper of OvchinnikovEl who argues that the decay of the SLS caused by these fluctuations diminishes with 
the increase of the mode amplitude. This statement, however, is based on a perturbational consideration 
which can not be applied rto the description of evolution of vibrations with a strong amplitude. 

Recently we showedli5li3 that the perturbational treatment of the effect of quantum and thermal 
fluctuations on a local mode associated with a defect atom in a lattice fails for the case when the 
amplitude of vibration is large. In fact, the two-phonon damping of the local mode, caused by cubic 
anharmonicity, behaves dramatically with the change of the amplitude: at definite "critical" amplitudes 
relaxation jumps take place being accompanied by a generation of phonon bursts. This effect indicates 
that the quantum and thermal fluctuations may dominate in the dynamics of strong vibrations with 
the energy of the mode being within a specific range. The strong field of the local vibratiop^causes 
the transformation of phonon operators and the increase of the number of phonons in timeU'113. This 
mechanism .of phonon generation by a local vibration is analogous to the mechanism of black hole 
radiatioiMM 

In this communication we extend the theoryEirliij to the SLS in a monatomic one-dimensional lattice 
(chain) with cubic and hard quartic anharmonicity. This case is of special interest since the monatomic 
chain with both anharmonicities is the simplest model for the investigation of the quantum relaxation of 
SLS in a perfect lattice. Hard quartic anharmonicity is a prerequisite of existence of the SLS, while the 
cubic anharmonicity stands for the two-phonon decay of the SLS. We note that SLSs in this model are 
stable in classical limit, since all harmonics of the mode are out of resonance with lattice phononsa. The 
odd SLS is examined and analytical nonperturbative solution of the problem and results of numerical 
calculations are presented. 

The potential energy operator in a monatomic one-dimensional lattice, which includes linear, and the 
first two nonlinear terms and takes account of the nearest neighbour interaction, has the following form: 

n r—2 

where U r is the operator of the longitudinal displacement of the n-th atom from its equilibrium position, 
K r are harmonic (r = 2) and anharmonic (cubic: r = 3, quartic: r = 4) springs. The operators U n 
satisfy the following equations of motion: 



d 2 U n 
dt 2 

r=2 



Y, [(Un+1 - Unf- 1 - (Un ~ Un-lY' 1 ] , (2) 
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K r = K r /M, M is the mass of an atom. We suppose that an SLS with the frequency uii < 2ujb is 
excited at the time t = at the sites n = and its neighbours (cod is the maximum harmonic frequency, 
lud = 2\/ K2). Anharmonic interactions are supposed to be weak, satisfying the condition <C Aq, 
where ~ ^h/ (2lue>M) is the amplitude of zero-point vibrations, Aq ~ K 2 /K4 is the amplitude 
of the self-localized vibration (from physical reasons the value of Aq should not exceed the value of the 
lattice constant d, i.e. K4 > K 2 /d 2 ). The condition a^ 1 <C Aq means that the characteristic energy of 
the SLS, being of the order of K2/K4 is much larger than the characteristic vibrational quantum hwo'- 

Kl/K 4 Kl' 2 MV* 
K= H^= 2HK, >>1 - (3) 

Note that the reversed dimensionless parameter /C _1 characterizes the degree of the quantum anhar- 
monicity of the lattice: it increases as the anharmonic term K4 increases and as the mass of the atoms 
M decreases. In quantum crystals He and Ne K, is of the order of 1, while in ordinary crystals (char- 
acterized by small amplitude of zero-point vibrations as compared to the lattice constant) it can reach 
many hundreds. In the problem under investigation K, plays an essential role determining the time-scale 
of the energy relaxation (see below) . 

In order to take account of the SLS decay we introduce the operators U n in the form 

U„ = u n (t) +£„ + q n , (4) 

where 

the "classical" displacements u n are supposed to be nearly periodic functions: u n (t) ~ A n {t) cos u>it, 
A n (t) and (,n(t) are, slowly changing with time, amplitudes and shifts satisfying the equations 

- uj?A n = K 2 (A n+1 - A n ) + 2K 3 (A n+1 ^ n+1 - !„£„) + 

+ 3 /? 4 (^±i + A n+I e n+1 - ^ - A^r) , (5) 

- A 2 A 2 



= k 2 n n+1 - e„) + k 3 (-^± + g +1 - -f - g 

+ K i {^^ + ^ +1 -^-^), (6) 

where we keep just the cos u>t term neglecting ~ cos Suirf , cos 5uit , . . .!~@ and introduce the new variables: 
A n = A n — An— 1, £ n — £ n — £ n -i- Variations of the A n and £ n in time, due to quantum fluctuations, 
are described by means of time-dependent operators q n (t), which satisfy the equations 

-Jp = W n+1 (t)(q n+1 - q n ) - W n (t)(q n - q n ^) . (7) 

Here W n (t) = K 2 + k n + 2w n cos w/t, where k n = 2K 3 £ n + 3K A (Al/2 + g*) and w„ = K 3 l„ + 3^ 4 A n f n 

determine the change of springs caused by the SLS. The terms ~ K 3 (q n — q„_i) 2 , K±A n (q n — g„_i) 2 , 
Ki^niQn ~~ <?n-i) 2 and K^{q n — q n -i) 3 are supposed to be small and they were neglected. It is important 
to take these terms into account when u>i > 2u>d, i-e. when the two-phonon decay under consideration 
is forbidden by the energy conservation law. 

The phonon hamiltonian is defined as H ph (t) = T(t) + W(t), where T = \ Y, n ^n and W = W0 + W1 + 
Wt are the operators of kinetic and potential energy, Wo = K 2 /2 ^2 n (q n — Qn-1) 2 , W\ — 1/2 J2 n ^n(<Zn — 
<Zn-i) 2 , Wt — cosuiitJ2 n w n {q n — q n -i) 2 ; K=l. Wi describes the stationary perturbation of the phonon 
subsystem by the SLS, while Wt takes account of the oscillatory time-dependence of the springs induced 
by the SLS. This is the-kipae-dependence of H (t) which causes the phonon emission by the SLS (phonons 
are generated by pairsti'EJ). 

Introducing properly chosen configurational coordinates X{ = ^2 n Si n q n and Q m = ^2 n s mn q n , one 

can diagonalize Wi and Wt- W\ = ^J^iViXf, Wt — cosojit J2 m v <nQ 2 n - Both W\ and Wt have one 
zero eigenvalue; all other rji are positive, while v m are sign-alternating, being symmetric with respect to 
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the sign change. The mode with zero eigenvalue is totally symmetric and therefore does not enter into 

The Hamiltonian H p h is diagonalized by introducing time-dependent operators as follows. First 
the time-independent Hamiltonian H — + Wq is diagonalized by applying the standard methods 



of local dynamicsE3Ej (6j$ is the Hamiltonian of the chain in harmonic approximation). Then the 

Hamiltonian H p h(t) is diagonalized by the method presented in0O'EI The timc-depcndcnt phonon 
Hamiltonian in the diagonal representation-allows one to find energy of phonons E p h (t) and the rate of 



= I(t). (8) 



phonon generation. The latter equals t( 

dEphif) = dEi{t) 
dt dt 



Here E l ~ (^M^- + J2 r K r^) is the energy of the SLS, 

r° dwY^ nLJtV " 1 Prn{u)Pm{ui ~ + 2n(u})) ^ 

Ju,-uJn ~f 8 |1 - V 2 G m (oj)\ 2 



I = 



denotes the intensity of the emission of phonons, p m = —$s(G m (uj)), the temperature factor is included 
by the multiplier (1 + 2n(u>)), G m (ui) = G m (w)G m (w — u>i), 

POO 

G m (w) = -i dte iut - et ([Q m (t),Q m (0)}), e^O, (10) 



— „iHtr> e ~iHt 



is the Green's function of the mode m, Q m (t) = e Q 

Use of the diagonal representation for the operator W\ allows one to easily express the Green's 
functions G m via the one-cite Green's function of the perfect latticeEil Gq°' ) = i/{uj^Juj 2 D — to 2 ) as follows 



n - M°) \ " VimOi'mSinSi'ni cos(n(n-n')) r 

" " h, (1- W 0T)(1-*G|!») 1 """ 



(0) 



(11) 



Here G| 0) = g[, 0) £ nn , S in S in , cos(fc(n - n')), a mi = Y Jn s mnS tn , cosk = 1 - 2u 2 /u 2 Dl uj 2 < lo 2 d . 

In the limit of small \v m \ one can neglect the term w^ 1 G rrl (w) in the denominatox-pf formula (0). In 
this limit (9) coincides with the corresponding formula of the perturbation theorjO. For realistic K 3 
and the considered values of loi between 1.3w£> and 2lud the term v^Gm is not small. At some uj\ = uj^ 
and uj both 1 — v 2 n ?R.(G m (oj)) and ^s(G m (ui)) turn to zero and the perturbation theory fails. Near such 
ujk I ~ \uji — uJk}^ 1 ~ \t — tk\ _1 / 2 , i.e a sharp burst of phonons is generated causing a relaxation jump 
(ijt corresponds to the time moment when uji — u>k)- 

Note that is not a linear function of the energy of the SLS. Therefore the relaxation of the SLS 
is always nonexponential, including the case when perturbation theory is applied. In the latter case 
Ei(t) ~ (ti —t) a , where a ~ 0.75 . . . 0.83, f; is the lifetime of the SLS which is finite in this approximation. 

As an example the quantum decay of an odd SLS is examined below. The properties of the SLS depend 
essentially on the dimensionless parameter S — \/K 2 j (K2K4). For realistic one-well pair-potentials this 
parameter has the value between 1 and a/32/9. For Ar-Ar and Kr-Kr potentials 5 = 1.37, for K-Br 
potential 6 — 1.31. For such values of 5 the odd SLS is well localized if u>i ~ I.Sojd- In this case SLS is 
localized on three central atoms: Ai = A_i ~ —Aq/2, \A n \ ~ 0.05|A |, n > 2; = — £_„ — £, n > 1, 
£0 = Olru. This allows one to take account of only Aq, Ai and £x = — £_i ~ £. The parameters Aq and 
£ of the classical SLS problem in this approximation are determined by the equations @ [n — 0) and 
© (n = 1): 



OI 

cjf = 3K 2 + 6K 3 £ + —K4A 2 , + 9K^ 2 , (12) 



= K 2 £ + K 3 (A 2 Q + e) + K 4 ^(y A 1 + • (I 3 ) 
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Introducing dimensionless amplitude a = y K^A^jKi and shift z — £y/ K4/K2, one gets 

« 2 = (14) 
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where z is determined by the real solutions of the equation 



3 81 z(5z 2 + fza+ j* 2 -l) 
' ~ 4 64 5Tfi (15) 



(here and below the units u)d = 1 are used). This is a cubic equation for z which for considered uji and 
8 has only one solution with a 2 > 0. This solution determines the parameters of the SLS in the classical 
limit. 

In the approximation considered, the SLS causes changes of only in central and next to central springs: 

fc = fc 1 = (16z<S+24z 2 + 27a 2 )/32, k 2 = k_ 1 =3a 2 /32 , 
wq = —w\ = a(S + 3z)/8 , W2 — —W-i — a8/8 . 

In this case W± and Wt can be diagonalized analytically. The eigenvalues T]i and v m and the components 
of the eigenvectors are the following: rji = 0, 



2 J V 2 

U3 = , «i,2,4,5 = ±^ (8^2 + 6w 2 + 4wi w 2 ± 
2(16w2 + 8w 2 w 2 + Wwlw! + $w\ + I2wlw 2 )^) * 
(H — h corresponds to v\, — h to u 2 , H — to V4 and to V5); 

Skn — — 7==== ) = Sk,l = 1 , S'fc2 = 1 — %/ ) 

S fc3 = -4 + 2 % //3 , S ki = 1 - , Sk5 = 1 , * = 2, 3 ; 
<Sfe3 = , SW = -1 + W/3 > -5fc5 = -1 , fc = 4, 5 . 
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Smfi ~ _ X * ~ 

Sinn / \~^ ~2 ' — / ^in : 

V ^mn n—2 

s i2 = l- Vi(ai - w 1 /w 2 + l)/wi , s i3 = 1 , s t4 = 1 + a j; , 
Si5 = 1 + Vi/w 2 , a { = Vi(vi +w 1 + 2w 2 )/(w 1 w 2 ) ■ 

Energy of the SLS equals 

Ei ~ ?iu D ]Csi; ei = Scofa 2 + z 2 {^ + | + (16) 

one sees that parameter /C determines the quantum scale of the SLS energy. Note that the intensity 
of phonon emission does not depend on K, (the unit of / in (^) is fvJ\^). Consequently, parameter K, 
determines the time-scale of energy relaxation. 

We performed calculations of the frequency dependence of the intensity of phonon emission by the SLS 
and of the time dependence of the intensity of phonon emission and the SLS energy (at the temperature 
T = 0). The calculation procedure was the following. First we fixed the initial SLS frequency (usually 
at oji = 1.98a;£)) and calculated z, a 2 , Ei, G(ui), I. The second step was to take a lower value of u>i and 
respectively calculate new values of parameters. We repeated this procedure until the frequency of the 
SLS got value 1.3. The time step was calculated using the relation At = —AEi/I. 
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In Fig. 1 the dependence of the intensity of phonon emission on the value of u>i is shown. Indeed, the 
intensity / has peaks ~ \uii — uJkl^ 1 at critical frequencies uik (this was confirmed numerically). This 
means that the SLS is unstable in the vicinity of ojk . The time dependences of the intensity of phonon 
emission and of the mode energy E\ are given in Fig. 2. 

In all cases the lowest considered frequency u>i — I.Slud corresponds to a small energy as compared 
to the initial energy for uii ~ 2u>r>. Consequently, the main part of the energy of the SLS with initial 
frequency lui ~ 2ojd is lost during the time considered (note that an extrapolation of Ei(t) to larger t 
also gives finite lifetime of SLS). Depending on the value of 5 and on the initial energy, this time can 
be either longer or shorter than that given by the perturbation theory. The relaxation law is essentially 
nonexponential: sharp bursts of phonons are generated at critical values of time tk when uii approaches 
uik causing relaxation jumps. The time dependence of the intensity I near the t k is I ~ \t-t h \~ 1 / 2 . The 
peaks mentioned are caused by the poles in the integrand function (^); they correspond to the emission 
of quasimonochromatic phonons and give sharp lines in the spectrum of generated phonons J(u>) (the 
integrand function in (||) summed over time) in Fig. 3. 

The SLS is not fully localized on three central atoms: it always has wings. Equations (4) and (5) give 
the following solutions for the amplitudes A n and the shifts £„ in the wings: 

A n ~ (-l)Moe-^H , ^~sign(n)^° e - 2 ^"l, n » 1 (17) 

where fj, = 2arch(w//wi)). It can be shown that taking into account these amplitudes and shifts leads to 

relaxation jumps of SLSs with oji close to 2u)d (2ujd — oji ^ 0.01). These relaxation jumps are weaker 
than those caused by the central atoms.. 

In the case 5 2 < 3/4 there may exisilj large-size SLS with wi close towi). The energy of such an SLS 
Ei ~ d>y/2{u)i — 1)K 2 /{K±{Z — 45 2 )) exceeds the unit (Tiujd) only if S 2 approaches 3/4 from below. The 
problem of quantum decay of this large size SLS is analogous to that of solitons in continuous media 
and it will be considered elsewhere. 

In conclusion, we examined the quantum relaxation of the odd SLS in a one-dimensional anharmonic 
monatomic lattice, caused by the zero-point fluctuations of the phonons and gave an analytical nonper- 
turbative solution of the problem. We found that the fluctuations can cause generation of very sharp 
bursts of quasimonochromatic phonons and very fast, step-wise jumps of the SLS energy. The rate of 
energy loss is on average of the order of the phonon quantum hujo per period of vibrations 2-k/ujd- The 
full relaxation time is determined by the quantum anharmonicity parameter /C; in chains with poten- 
tials corresponding to the quantum crystals, this time is of the order of a period of vibrations, while in 
lattices with potentials of ordinary crystals it may reach thousands of periods. 
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Figure Captions 

Figure 1. Frequency dependence of the intensity of emission of phonons by the odd SLS.Thc peaks corre- 
spond to the regions where the SLS is unstable withrespect to quantum fluctuations; 5 = {K%/ K 2 K i ) 1 / 2 



Figure 2. Time dependence of the SLS energy E\ (thick dashed line) and intensity of phonon emission 
by the SLS / (thick solid line) ; decay of the SLS energy calculated with perturbation theory (thin solid 
line). Sharp peaks in the I(t) dependence describe emission of phonon bursts by the SLS. E\ = EiK^/K^, 
K, = K\j{K^hiOD)- I is measured in units of Huj%. 

Figure 3. Spectrum of generated phonons J(u>) (in a.u.). 
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